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The paper is about de-noising procedures aimed on events recognition in signals from a distributed 
fiber-optic vibration sensor system based on the phase-sensitive optical time-domain reflectometry. 

We report experimental results on recognition of several classes of events in a seismic background. 

A de-noising procedure uses the framework of the time-series analysis and Kolmogorov-Zurbenko 
filtering. We demonstrate that this approach allows revealing signatures of several classes of events. 


I. INTRODUCTION 

Real-time monitoring systems with the use of dis¬ 
tributed fiber-optic vibration sensor systems based on the 
phase-sensitive optical time-domain reflectometry tech¬ 
nique [1-6]. have a fascinating prospective for appli¬ 
cations. Examples include control for access on pro¬ 
tected areas {e.g., securing national borders), oil and gas 
pipelines, communications lines, and structural health 
monitoring [7-9]. 

The core of such system is phase-sensitive optical time- 
domain reflectometry technique, which has sufficiently 
high sensitivity and spatial resolution [10-18]. A main 
feature of this type of reflectometry is a sufficiently large 
coherence length of the employed optical pulse. Signals 
reflected from centers of the Rayleigh backscattering ex¬ 
hibit the coherent summation of their complex wave am¬ 
plitudes. 

On the one hand, monitoring systems based on the 
phase-sensitive optical time-domain reflectometry are 
sensitive enough to register sufficiently small fluctuations 
[6]. On the other hand, this means that an algorithm, 
which allows revealing the nature of fluctuations, should 
supplement such systems. Indeed, the crucial problem 
here is to reveal: are these fluctuations caused by natu¬ 
ral changes of background or by any kind of activates? 
In other words, a highly non-trivial problem of de-noising 
comes to the fore. 

To make a decision about the nature of fluctuations 
one can continuously analyze signals from the system in 
time or frequency domains. Due to a sufficiently complex 
structure of signals from the monitoring system this prob¬ 
lem is rather challenging. It has been extensively studied 
during last decade (see [18-20] and reference therein). 
However, at this moment there is no universal solution 
for events recognition problem for vibration sensor sys¬ 
tems based on the phase-sensitive optical time-domain 
reflectometry technique. 

Mathematically speaking, signals from distributed 
fiber-optic vibration sensor systems are time series [21]. 
Then for their analysis various types of time series anal¬ 
ysis can be applied. In this paper, we present exper¬ 
imental results on an application of the Kolmogorov- 
Zurbenko filtering [22-26] for signal de-noising in a flber- 



FIG. 1. Fiber-optic distributed sensor system based on the 
phase-sensitive optical time-domain reflectometry: setup for 
collecting of experimental data. The sensing element of the 
system is the standard single mode fiber (fiber-optic cable). 

optic distributed vibration sensor system based on the 
phase-sensitive optical time-domain reflectometry. In the 
considered case, the main goal of the system is a con¬ 
trol for access on protected areas. Consequently, we are 
confronted with the problem of events recognition {e.g., 
human passage, human group passage, or car travel) in 
a seismic background. 

The paper is organized as follows. In Sec. II, we de¬ 
scribe our setup for collecting experimental data and pa¬ 
rameters of the fiber-optic distributed vibration sensor 
system. In Sec. Ill, we describe the basic de-noising 
procedure based on the Kolmogorov-Zurbenko filtering. 
In Sec. IV, we present experimental results on applica¬ 
tion of the Kolmogorov-Zurbenko de-noising procedure 
to measured signals. In Sec. V, we give our conclusion. 


II. SETUP FOR COLLECTING 
EXPERIMENTAL DATA 

The fiber-optic distributed vibration sensor system for 
collecting of data is located on the Bauman Moscow State 
Technical University polygon in Moscow Region. The 
setup is presented on Fig. I: I is the primary light source 
(laser), 2 is the acousto-optic modulator, 3 is the circu¬ 
lator, 4 is the fiber-optic sensor, 5 is the preamplifier. 
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6 is the optical filter, 7 is the detector, 8 is the ADC 
converter, 9 is the programmable logic device, 10 is the 
computer. 

Probe signal has the wavelength 1550 nm, probe pulse 
of duration 200 ns, and the signal from the semiconductor 
laser of power 300-500 mW is launched into the standard 
optical fiber. Probe signal has ultra-narrow line width, 
which is less than 1 MHz. In our experiments, length of 
the optical fiber cable 1 is approximately 50 km. 

The idea of its work can be presented as follows. In 
case of a vibration impact, the intensity of backscattering 
light changes according to the level of the impact. Cir¬ 
culator is used for launching of the probe signal into the 
optical fiber and the backscattering signal to a detector. 
Signals in the system are sum of all scattered signals dur¬ 
ing time of pulse with taking into account their phases. 


III. SIGNALS FROM THE SYSTEM: EVENTS 
RECOGNITION 

A typical result of the measurement using the setup 
is presented in Fig. 2a. This figure presents an overlap 
between signals measured by our setup (Sec. II) on the 
region of the cable with length 0.5 km during one second. 

As it was mentioned above, the crucial challenge is to 
recognize any kind of deliberate activity in these signals. 
It is clear that this problem can be essentially divided 
on two related sub-problems. The first part is to regis¬ 
ter the event in a (seismic) background. It is seen from 
Fig. 2b that an event can be registered by the system via 
measurement of the difference between signals in neigh¬ 
boring moments of time. From such a procedure, one 
can find, e.g., position of the event in the cable. How¬ 
ever, due to natural fluctuation of the background such 
simple procedure leads to a sufficiently large number of 
false positives. 
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FIG. 2. An overlap of typical signals measured by the system 
on the region of the cable with length 0.5 km during one 
second (left) and the maximum deviation (right): a signature 
of an event. 

The second problem is to classify registered event. In 
this work, we are guided by the following basic classifi¬ 
cation of events: (i) single event, which is localized in 
space and in time; (ii) single event, which is delocalized 
in space and localized in time; (iii) single event, which 
is localized in space and delocalized in time; (iv) single 
event, which is delocalized in space and delocalized in 
time. 

For our setup an important part is using of preliminary 
tests. These test consist of collecting of a large number 




of experimental data corresponding to background and 
typical types of activates. Preliminary test have been 
organized on the polygon at similar experimental condi¬ 
tions (importantly, the same climate conditions). From 
preliminary tests, we use a sufficiently large number of 
experimentally measured signals of type (i)-(iv) to obtain 
minimum, maximum, and average characteristics values 
(characteristic time scales of events and characteristic 
length scales) as well as parameters of the background. 



FIG. 3. Measured signals [27]: (a) background signals exper¬ 
imentally measured by the system; (b) sequence of space local¬ 
ized short-term events with additive noise form background; 
(c) sequence of space delocalized short-term events with ad¬ 
ditive noise form background; space localized long-term event 
with additive noise form background; (e) space delocalized 
long-term event with additive noise form background; results 
of the de-noising procedure (f)-(j) for the signals (a)-(e). 


The background noise is the waterfall form, z.e., inten¬ 
sity as a function of time (t) and position in the fiber¬ 
optic cable (1), is presented in Fig. 3a. A sequence of 
space-time localized events is presented in Fig. 3b. Fig. 
3c shows a sequence of delocalized in space and localized 
in time events. In Fig. 3d, one can see a sequence of 
localized in space and delocalized in time events. Fig. 3e 
shows a sequence of space-time delocalized events. 
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IV. BASIC ALGORITHM 

Signals from a distributed fiber-optic vibration sensor 
systems are time series. Therefore, for their analysis var¬ 
ious types of time series analysis can be employed. One 
of the possible solutions is to use statistical apparatus of 
the time series analysis. 

For instants, an event can be detected in the back¬ 
ground in the difference in the value of signals between 
time moments t = j and t = j + 1. Toward this end, 
we suggest to implement the spatial low-pass filter of the 
Kolmogorov-Zurbenko type [22-26]. Similar approaches 
to de-noising procedure based on moving average calcu¬ 
lation have been recently discussed [28, 29]. 



FIG. 4. 3D visualization of the signal measured form the 
system after the de-noising procedure based on the basic al¬ 
gorithm (l)-(4) and classified by the system as sequence of 
space and time localized events, which can correspond to hu¬ 
man passage along the cable. 

In general, the Kolmogorov-Zurbenko filter is a series 
of iterations of the moving average filter. The first iter¬ 
ation of Kolmogorov-Zurbenko filter is an application of 
the moving average filter over a given process. 

(m-l)/2 

KZ[x(i,0L,,=i= E + W 

s= —(m—1)/2 

Here, X{l^t) is the signal (real-valued time series), m 
is the time window, and k is the filter order [22]. The 
iteration process of a simple operation of moving aver¬ 
age is very computationally efficient. The Kolmogorov- 
Zurbenko filter has been employed in investigations of cli¬ 
mate fluctuations and seminal studies of the turbulence 
problem [22-26]. 

Using the de-noising procedure based on the applica¬ 
tion of the Kolmogorov-Zurbenko filter (1) to signals, the 
sensor system can detect an event in the background. In¬ 
deed, if the moving average deviation of the signal A(/, t) 

= ( 2 ) 


exceeds the critical value, then a potential event is de¬ 
tected. The key question for implementation of such an 
approach is about determination of the critical value (2). 

We suggest use an adaptive critical value calculation 
in a rather straightforward way. Let us consider the def¬ 
erence between two signals with a time window j in the 
following form: 

5i{M) = \X{t + j,l)-X{t,l)\. (3) 

By integrating this difference over the cable length 

JSj{Al)dl = S* (4) 

one can obtain the critical value, which characterizes the 
integrated difference between two signals from the sys¬ 
tem in neighboring countdowns of time. On practice, we 
obtain critical values on the bases of preliminary tests of 
the system (see Sec. III). 

This critical value is continuously measured by the 
monitoring system during preliminary tests. Finally, 
from the de-noising procedure with taking into account 
the critical value, we obtain the de-noised signal as fol¬ 
lows: 



KZ[X{t,l)], 

0, otherwise. 


( 5 ) 


We optimize the de-noising procedure with respect to 
parameters of the Kolmogorov-Zurbenko filter: size of 
the time window and order of the filter. This allows 
obtaining the maximal level of recognized events. 

The results of the application of the de-noising pro¬ 
cedure (l)-(4) for experimentally measured background 
signals with embedded events [27] of type (i)-(iv) are pre¬ 
sented in Fig. 3. In spite of simplicity of used de-noising 
procedure, events can be clearly detected. Furthermore, 
patterns of different classes of events can be reliably re¬ 
vealed (See Fig. 4). 


V. RESULTS AND CONCLUSION 

Controllability of a workflow is the most important re¬ 
quirement for its efficient implementation. Then a prob¬ 
lem of a design of real-time monitoring systems for non¬ 
destructive testing has a paramount importance for tech¬ 
nological processes in a broad class of applications in in¬ 
dustry. Mainly, the goal of a monitoring system is to 
measure and analyze specific characteristics of a process, 
which allows to make a decision about the correctness of 
its functioning. 

In this work, we have investigated the basic algorithmic 
solution for recognition of events in seismic backgrounds. 
Performance of the suggested algorithmic solution has 
been demonstrated on experimentally measured signals. 
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